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Abstract 



We considered Q-state Potts model on Bethe lattice in presence of external magnetic 
field for Q < 2 by means of recursion relation technique. This allows to study the phase 
transition mechanism in terms of the obtained one dimensional rational mapping. The 
convergence of Feigenabaum a and 6 exponents for the aforementioned mapping is in- 

■ vestigated for the period doubling and three cyclic window. We regarded the Lyapunov 
^ I exponent as an order parameter for the characterization of the model and discussed its 

dependence on temperature and magnetic field. Arnold tongues analogs with winding 
numbers w = 1/2, w = 2/4 and w = 1/3 (in the three cyclic window) are constructed for 
, Q < 2. The critical temperatures of the model are discussed and their dependence on Q is 

I investigated. We also proposed an approximate method for constructing Arnold tongues 

via Feigenbaum 6 exponent. 

O ' 
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1 Introduction 

The integer Q-state Potts model is a gener- 
alization of the Ising model to more-than-two 
components. Some particular cases of the model 
(different integer values of Q) played an es- 
sential role in the theory of phase transitions 
and critical phenomena [HIS]- For Q — )■ 1 it is 
associated with the percolation (connectivity) 
model [2111]. Both dilute [5] and "pure" Potts 
models can be also represented by Kasteleyn- 
Fortuin (KF) clusters fB]. The partition func- 
tion of the Potts model is well defined for any 
non-integer Q > |7j. Some interesting special 
cases arise in the limit Q — )■ 0: e.g. the limit 
Q]v = e^'^ — 1 — )■ with u = ^ held fixed gives 
rise to a model of spanning forest [HIIS] . 

Potts model can be also realized in exper- 
iments. In particular, Q-state Potts lattice 
model can be used for description of cold de- 
naturation of proteins in solvents p^fTT]. Re- 
search of experimental realizations of the model 
shows that in the case < Q < 1 it is in the 
same universality class of the transition in gela- 
tion and vulcanization processes in branched 
polymers [12]. Besides, many physical pro- 
cesses like the resistor network, dilute spin glass, 
self organizing critical systems can be formu- 
lated in terms of the model, when Q < 2 [13^ 

[ig. 

Since the Potts model generally is not ex- 
actly solved (there isn't any solution of the 
model in dimensions > 2 in non-zero mag- 
netic field), approximate methods are required. 
An important advance was made with the in- 
vention of the Swendsen-Wang (SW) and 
Chayes-Machta (CM) [16] algorithms for sim- 
ulating the ferromagnetic Potts model at pos- 
itive integer Q and the random-cluster model 
with any real Q > 1 respectively (for a sum- 
mary see [17]). No analytical results can be 
obtained in these cases. 

In this paper we used the Bethe approx- 
imation. One of the advantages of this kind 
approach is that some models can be solved 
here exactly with the help of the recurrence 
relation technique [18] and analytical results 
can be obtained. The method is more reli- 



able than mean-field calculations [19]. This 
approach can be also applied to the general- 
ized Bethe lattice (known as Husimi lattice), 
which can be used for investigation of models 
with multisite interactions [20] . It was used for 
the study of magnetic properties of the solid 
^He films ^IJ. In the case of recurrence re- 
lation technique, the properties of the model 
are associated with the behavior of the iter- 
ation sequence {xn} of a nonlinear mapping, 
which mostly exhibits period doubling cascade, 
chaos, p-cyclic windows. Here the Lyapunov 
exponent serves not only a good order param- 
eter, but also gives relevant information about 
the geometrical and dynamical properties of 
the model's attractors [22]. The phase struc- 
ture of the model can be investigated with the 
help of mapping's Arnold tongues [25] . 

The concepts Scaling and Universality have 
played an essential role in the description of 
statistical systems [21]. Particularly the be- 
havior of famous Feigenbaum exponents [251 
[2B] for one dimensional rational mapping de- 
scribing statistical model on Husimi lattice was 
investigated [27] . 

This paper is organized as follows. In the 
next section we present the recursion relation 
and phase structure of the Potts model on the 
Bethe lattice. In Sec. [3] the universality of 
Feigenbaum exponents for the obtained one- 
dimensional rational mapping is derived (for 
Q < 2). The dependence of the Lyapunov 
exponent on external magnetic field and tem- 
perature is given in Sec. HI In Sec. E] Arnold 
tongues analogs are constructed for Q < 2. In 
Sec. [6] we investigated the three periodic win- 
dow of the mentioned above mapping. The 
critical temperatures and an approximate me- 
thod for constructing Arnold tongues through 
Feigenbaum 6 exponent are introduced in Sec. [71 
Finally we summarize our results in Sec. El 
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2 Rational Mapping for the 
Potts Model and Phase 
Structure 

The Q-state Potts model in the case of two-site 
interactions in presence of external magnetic 
field is defined by the Hamiltonian 

n = -jJ2S{a„a,) - if^5(a„Q). (1) 

The first sum in Eq. ([1]) goes over all nearest- 
neighbor pairs and the second one over all sites 
of the lattice (J < is the anti- ferromagnetic 
and J > the ferromagnetic case). The par- 
tition function and single site magnetization is 
given by 



Z 



(2) 



M = {Siao,Q)) = Z-'J2^(''o,Q)e 



(3) 



where fcs is the Boltzmann constant (we will 
set fee = 1). The Bethe lattice can be sepa- 
rated into 7 identical branches by cutting apart 
at the central point. The partition function 
can be written as 



Zn 



exp{^-5iao,Q)}[9n{cTo)V, (4) 



where ctq is the central spin and (7„(c'"o) is the 
contribution of each lattice branch. Following 
well-known procedure [T8l[20l|21] , we obtain 



fix) 



where 



+ (e^ + Q - 2)x^-^ 
e^ + {Q-l)x^-^ 

9n{(T ^ Q) 



(5) 



(6) 



gn{cT = QY 

Equation (E]) is called Potts-Bethe mapping. 
Knowing recursion relation we can calculate 
the magnetization of the central spin: 



M„ = (5(cro,g)) 



H_ 

e T 



(7) 



In real statistical systems the bifurcation 
points of mappings hke Eq. ([5]) correspond to 
phase transition points of second order. For 
systems with Q < 2 and anti-ferromagnetic 
coupling (J < 0), the dependence of M mag- 
netization on the field H is complicated: full 
range of period doubling, chaos, p-cyclic win- 
dows can be observed here. Figure [T] shows 
bifurcation diagrams for different values of pa- 
rameters. 

In terms of the mapping in Eq. (|5]) the 
area, where M is single-valued function of H 
[AB and CD in Fig. HJa), AB in Fig. dJ^b)], 
the recursion sequence {x„} of the mapping 
converges to one stable fixed point xq (phase 
without sub-lattice structure). The area af- 
ter first bifurcation point [BC in Fig. [D^a), 
C1BC2 in Fig. dl^b)], sequence converges 
to two stable points. Therefore here we have 
two values of magnetization, which should be 
explained, as an arising of two sub-lattices of 
anti-ferromagnetic order. The areas between 
consequent bifurcation points are described by 
a sequence of n fixed points (existence of n sub- 
lattices), corresponding to various modulated 
phases with finite period (commensurate mod- 
ulated phases), which also appear in p-cyclic 
windows. 

3 Universality of 

Feigenbaum Exponents 
for Potts-Bethe Mapping 

Let Tn be the value of the parameter r of a map- 
ping at the period doubling bifurcation point, 
Too the value of r from which the chaotic be- 
havior ensues. It turns out that the values of 
r„ satisfy the following scaling: 

''^n = ^oo — const6~^, n ^ 1. (8) 



If dn is the distance between the point x* for 
which fr{x) is extremal and the nearest point 
on the 2" cycle (Fig. [2]), then 

= -a, n > 1. (9) 



d. 



n+l 
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(a) (b) 




Figure 1: The plot of M (magnetization) versus external magnetic field H for 7 = 3, (a) Q = 2.5, 
J = —1, T = 0.5 (the areas AB and CD correspond to phase without sub-lattice structure, BC to 
antiferromagnetic phase, B and C are the points at which the phase transition between these phases 
takes place); (b) Q = 1.1, J = —1, T = 0.5 {AB correspond to the magnetization of the paramagnteic 
phase, BCi and BC2 to magnetization of two sub-lattices in the antiferromagnetic phase, B is the 
point of phase transition between paramagnetic and antiferromagnetic phases, C\ and C2 the points 
of transition between antiferromagnetic and four-periodic modulated phases). 



The quantities a in Eq. ^ and 5 in Eq. 
are Feigenbaum exponents 



5 = 4,6692016091... 
a = 2,50290787050... 



(10a) 
(10b) 



If Rn is the value of the parameter at which the 
line X = X* intersects 2" periodic cycle (Fig. [2]), 
then 

Rn = Roc- const' d-"". (11) 

The exponents 6 and a are universal, i.e. 
Eqs. (IE]), (E]) and ( fTTj) are true for wide variety 
of mappings, a and 6 having the same values 
as in Eqs. (ITnall and (flObl) [25]. For the const 
and const', presented in Eqs. ([8]) and (ITT!) , they 
depend on family of reflection functions. 
Introducing 5„ 



R„ — R 



n—l 



taking into account that 

S = lim 5, 



Rr. 



(12) 



(13) 



and the fact that 



■iRn ' 



X*) 



X 



{f^"\x) is the n — th iteration of f{x), 



I.e. 



etc.), we can calculate Feigenbaum exponents 
for a mapping f{x). In our case 



7-1 



CT + (cT + Q — 2)x' 
eT' -67 + {Q - l)x'5'-i 
r + (e^ + g - 2)xT-i 
r • + (Q - l)xT-i ' 



where 



R2 h Rj 



(15) 



(16) 




Figure 2: The period doubling process (schemat- 
(14) ically). 



Mapping f{x) is extremal for x* = 0. 



f('\x) = f[f{x)], /(2^)(x) = /W(a;) = /(2)[/(2)(a;)], 
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The values of i?„ and r„ must satisfy fol- 
lowing condition (see Fig. [2]): 



n < Ri < r2 < R2 < r3 < ... 



(17) 



We start solving Eq. iHM from n = 0, putting 
values of i?„ in order (JTTll . Besides, one can 
see, that 

lim R^ = R^ = Too- (18) 

n—>-oc 

Our numerical calculations for J = —1, 
T = 1,7 = 3, Q = 1.1 and Q = 0.8 are shown 
in Table [H 

We find that values of a and 6 converge to 
the Feigenbaum exponents for the Potts-Bethe 
mapping, which describes physical properties 
of a real statistical systems. The convergence 
in the case < Q < 1 is faster than in the case 
1 < g < 2 (see Table [I]). 

For const', presented in Eq. ffTTj) . ini the 
case Q = 0.8 and Q = 1.1 we obtained the fol- 
lowing values: const' = —2.682... and const' = 
—5.034... respectively. 

4 Lyapunov Exponents 
for One Dimensional 
Rational Mapping 

It is well known that under action of the map- 
ping Xn+i = f{xn) two nearby points can be 
dispersed. The Lyapunov exponent \{x) char- 
acterizes the degree of the exponential diver- 
gence of two adjacent points. The exact for- 
mula for the Lyapunov exponent is: 



A(x) = lim lim — In 

71— >oo e— >-0 n 



/(")(x + e) - /(")(x) 



examine the behavior of the magnetization at 
fixed external magnetic field, when tempera- 
ture is varied (see also Sees. [S]and[7]). 

Figures |3]^a) and (b) indicates, that at fixed 
H the chaotic region is richer than at fixed T: 
a large variety of different p-cyclic windows can 
be observed in Fig. MJo)- 

One of the major properties of the Lya- 
punov exponent is that it is equal to zero at 
the bifurcation point. This property can be 
easily obtained from the fact that the bifurca- 
tion point corresponds to a neutral point of the 
mapping. 

5 Arnold Tongues 

Values of parameters (external magnetic field 
and temperature) at phase transition (bifurca- 
tion) points can be found as a neutral point of 
the mapping: 



fix) 



(20) 



lim — In 

n— >-oo n 



dx 



(19) 



We regard the dependence of A(x) for the 
Potts-Bethe mapping [Eq. ([5])] on the external 
magnetic field H, fixing Q, temperature T and 
strength of interaction J [Fig. EJ^a)]. 

One can also study the dependence of the 
Lyapunov exponent on the temperature T, fix- 
ing H, Q, and J [Fig. Et^b)]. This will help to 



Different values of ip correspond to different 
types of bifurcation: 

1. (f = 2iTn. This case corresponds to saddle- 
node or tangent bifurcation (type I intermit- 
tency) [28l[29]. It can be observed e.g. in the 
case of the logistic map as ]?-cyclic window in 
the chaotic regime (see also Sec. [6]). 

2. ip = TT + 27m. Here we have bifurcation 
corresponding to period doubling [3Dj (type III 
intermittency) . 



3. A pair of conjugate complex eigenvalues of 
the Jacobian in the case of multidimensional 
mapping corresponds to the Hopf-bifurcation 
(type II intermittency) which introduces new 
basic frequencies in the system at the bifurca- 
tion points [3T] . 

At the first bifurcation point in period dou- 
bling regime we have: 



f{x) = X 
fix) = -1. 



(21) 
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(a) (b) 




Figure 3: The plot of the Lyapunov exponent for the Potts-Bethe mapping [Eq. dS])] (a) versus 
magnetic field H ior Q = 0.8, J = —1, 7 = 3, T = 2; (b) versus temperature T for Q = 1.2, J = — 1, 
7 = 3, F = 0.2. 



Eliminating x from Eq. ( 12T|) for the mapping 
in Eq. (jS]) and solving obtained equation in 
order to H, we obtain two branches: 



/ 



H = T-\n 



1 



.8(-l + Q) 

oJ/T 



_ 3J 

e T 



3J 

-6e T 



X (g - 2) - e^/^ (6Q + M - 6)(Q - 2) 

4 J 

-3e^ + (Q- 1)(Q + M- 1) 
-e^(3(Q-2)Q + M + 6 

and 
i7 = T ■ In 



(22) 



3J 

-e ' T 



3J 

6e T 



+ 

X (g - 2) - e'^/^(6Q - M - 6)(g 

4 J 

-3e- + (g-i)(g-M-i) 

+e^(-3(g-2)g + M-6 
where 



2) 



(23) 



X 



geJ/TQ _ g _ I8e^/^ + 9e¥ + l. 



(24) 



Let us consider two cases: 

1. g > 1. In this particular case two branches 
in Eqs. ( 122|) and (123!) define the hne separat- 
ing paramagnetic (P) and anti-ferromagnetic 
(2AF1, i.e. 2^ = 2 periodic) phases. 

2. g < 1. This case is more complicated. 
The reason is following (here and further 7 = 
3): when g < 1, Potts-Bethe mapping has 
singular points (see Fig. H]). 
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Figure 4: The plot of Potts-Bethe mapping f{x) 
[Eq. ©] for the case Q = 0.8, J = -1, 7 = 3, 
T = 2,H = 2. 

As one can see from Fig. HI the behavior of 
the mapping becomes sensible to initial point 
of the iteration: stable point can fall into area 
I or II . 

But it turns out, that only I has phys- 
ical meaning: in case g < 1, magnetization 
can be negative, which has not physical inter- 
pretation [see Eq. Since x„ converges to 
stable Xq point, therefore the criterion M being 
positive is 



H_ 

e T 



(25) 



1-g 

One can easily show that in area II condi 



tion is violated, hence here M < 0. So, 
in our further investigations we will assume, 
that stable points of the mapping are in area 
I. This case corresponds to the Eq. ( 123|) . It 
is also of interest to find the line separating 



Figure 5: Arnold tongues with winding numbers w = 1/2 and w = 2/4 for Potts model for Q = 1.2, 
J = — 1, 7 = 3 (the insert shows details in the area H > Q). Point C describes the upper bound of 
the temperature (Tc = 5.260340). 



anti-ferromagnetic (2AF1) and four - periodic 
(2M2, i.e. 2^ = 4 periodic) modulated phases. 
This means to find the values of T and if, 
where the second bifurcation occurs. The pro- 
cedure is the same with the only difference be- 
ing that here second iteration of the mapping 
(/(^)(x) = f[f{x)]) looses its stability: 

{f'"^^^)=^ (26) 

The area between two curves, obtained from 
Eqs. f l2T|) and fl26|) (anti-ferromagnetic phase), 
is the Arnold tongue analog with winding num- 
ber w = 1/2. Using the same technique we can 
also find the line, separating four - periodic 
(2M2) and eight - periodic (2M3, i.e. 2^ = 8 
periodic) modulated phases, i.e. the values of 
T and H, where the third bifurcation occurs. 
The condition will be: 

1 (/(^)(a:))' = -l. ^2^) 

The between curves defined by Eqs. fl26l) 
and f l27|) is the Arnold tongue analog with wind- 
ing number w = 2/4 (the area of existence 
four - periodic modulated phase). The result 
is given in Figs. |5] and [61 

It was already mentioned in Sec. [2] that 
convergence of Feigenbaum exponents in areas 



< Q < 1 and 1 < Q < 2 are different. As 
one can see from Figs. |5]and[6l Arnold tongues 
are different in these areas, too. 



a 





/W=l/2 


W=2/4 





1 2 3 4 5 



Figure 6: Arnold tongues with winding numbers 
w = 1/2 and w = 2/4 for Potts model for Q = 0.8, 
J = -1,7 = 3. 



M magnetization can be regarded as a func- 
tion of T, fixing the value of external magnetic 
field. Fig. |5] indicates, that lowering H first 
two periodic cycle appears, corresponding to 
P-2AF1 transition. When the line H inter- 
sects the upper curve of the Arnold tongue 
with w = 2/4, one finds a bubble on each 
branch [Fig. Cl^a)]. Here we have P-2AF1-2M2- 
2AF1 transitions. After the line H intersects 
the lower curve, we have another two bubbles 
which points on P-2AF1-2M2-2M3-2M2-2AF1 
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Figure 7: The plot of M (magnetization) versus temperature T for Q = 1.2, J = —1, 7 = 3 (a) 
H = 1.3; (b) H = 0.2. 



transitions. With further reduction of external 
magnetic field new bubbles are formed as parts 
of the old ones and for still lower if's we reach 
a region of chaotic behavior [Fig. El^b)]. 

6 Three Periodic Window 

In this section we will investigate the three 
periodic window of the Potts-Bethe mapping. 
Some definite values of parameters T and H 
form a line in the parameter space, which sep- 
arates the chaotic and three periodic regimes. 



(3) 



\ 1 \'-^ 






/ r" 






/ 1 "i^"" 






-4 / 1 1 / 
/ 1 1 /^-^ 

/ \l/ -LO 

/ -1.5 
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Figure 8: The Plot of /'^^•'(x) [third iteration of 
the Potts-Bethe mapping in Eq. ([5])] for Q = 1.1, 
J = — l,7 = 3at one point of the saddle-node 
bifurcation line defined by Eq. ([28]): T = 1.5915, 
H = -1. 

Figure IS presents the third iteration of the map- 
ping f{x) [Eq. ([5])] for one point from that line 
(as one can see the mapping has three fixed 
points here). On this line the saddle-node bi- 
furcation takes place, i.e. the aforementioned 



line can be found from the following condition 
(see Sec. [5]): 



/(3)(a;) = X 

(/(3)(a;))' = 



1. 



(2^ 



Consequent bifurcations correspond to pe- 
riod doubling, which leads to appearance of 
stable 3*2" periodic cycles. Following the tech- 
nique described in Sec. [5l we can find the line 
separating three and six periodic cycles (three 
and six periodic modulated phases) from the 
following system of equations: 



/(3)(a;) = X 

(/(3)(x))' = 



-1. 



(29) 




Figure 9: Arnold tongue with winding number 
w = 1/3 for Potts model for Q = 1.1, J = — 1, 
7 = 3 (the insert shows details in the area H > 0). 
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(a) (b) 




Figure 10: The plot of M (magnetization) versus temperature T in the three periodic window for 
Q = 1.1, J = -1, 7 = 3 (a) F = 1.3; (b) H = 1.2. 



Curves found from fl28l) and fl29l) form ArnoL 
tongue with winding number w = 1/3 (Fig. |9]). 
It corresponds to the area of existence of the 
three - periodic modulated phase (3M0, i.e. 
3 * 2° = 3 periodic). 

One can see from Fig. [9] that in the area 
H > the mapping presents interesting be- 
havior. Firstly, when the line H intersects 
only the upper curve of the Arnold tongue, two 
edges of the window are plainly distinguish- 
able: the saddle-node bifurcation takes place 
on both edges and the window is presented 
with only 3-periodic cycle [Fig. [TOTa)]. 

Transition between chaotic state and three 
- periodic modulated phase (3M0) occurs. Sec- 
ondly when the line H intersects lower line of 
the Arnold tongue, 3*2^ = 6 periodic cycle ap- 
pears in a form of a bubble, which corresponds 
to six - periodic modulated phase (3M1, i.e. 
3 * 2^ = 6 periodic). This indicates 3M0-3M1 
transition. It is obvious that if we continue 
lowering the values of H new bubbles will ap- 
pear and ultimately the chaotic region in the 
window will be reached [Fig. [TUT b)]. However, 
for H < the saddle-node bifurcation occurs 
only at one edge of the window. For thorough 
illustration of the effect, we also present the 
plots of Lyapunov A exponents in Fig. [11] 

We also investigated the behavior of Feigen- 
baum a and 6 exponents for three periodic win- 
dow. In comparison with the period doubling, 
the values of i?„ will be found from the follow- 



ing condition: 

fSf\^*)=^*- (30) 

Here again a and 6 converges to the values 
in Eqs. fllOal) and (jlObp . For example in the 
case J = —1, T = 1, 7 = 3 and Q = 0.8 for n = 
10, 6 = 4.669160924 and a = 2.502899839. 
For the const' in Eq. (iTTl) for the same values 
of parameters we obtained the following value: 
const' = —0.0167... These results confirm once 
more the universality hypothesis for rational 
mapping, which describes Q-state Potts model 
on Bethe lattice. 

7 Arnold Tongues 

through Feigenbaum 
Exponents and Critical 
Temperatures 

Constructing Arnold tongues for rational map- 
pings like Eq. ([5]), which describe real statisti- 
cal systems, is a very complicated and labori- 
ous procedure. But it turns out that universal- 
ity hypothesis (see Sec. [3]) is in deep relation 
with this problem. 



It can help in finding a simpler way for con- 
structing Arnold 

tongues, knowing only one line in parameter 
space, on which some kth bifurcation takes 
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place. Here will be the value of parameter 
[see Eq. f lTSj) ] at mth bifurcation point. From 
the scaling relation for and Eq. f lTSj) we can 
obtain in terms of {k < m), i?oo and 5: 



while Roo from Eq. (flTjl : 



5-1 



(31) 



(32) 



The convergence of the Roo to its real value 
is dependent on the speed of convergence of 6 
exponents (for Q>l?T, = 4or5 and in the 
area Q<ln = 5or6is enough). 

Since ri can be found analytically in con- 
trast with bigger n's, it would be convenient 
to take k = 1, and m > 1 [we already have ri, 
after solving Eq. (!2T1)]: 



ri + RooS - Roo J. n foo \ 

r2 = , for m = 2, (33aj 



5 

RqqS^ + n-Ro 
5^ 



for m = 3. (33b) 



For obtaining Arnold tongue analog with 
winding numbers w = 1/2 and w = 2/4 we 
have to find r2 and r^ for different T's and 
if's. The problem is that Eq. ([8]) is true for 
big n's. However, Eq. ([8]) remains true if we 



replace const with c„. In this case 



r2 



Cl 



+ Roc5 



2 _|_ C3(ri-_Roo) 



for m 



(34a) 



RooS 



Cl 



52 



for m = 3. (34b) 



As one can see, the value of r2 and rs from 



Eqs. dMID and fl33bl) closer to its real 



value in Eqs. fl34al) and (]34b|) as closer ^ and ^ 
is to 1. Values of ^ e.g. for investigated in the 
previous section Q = 0.8 and Q = 1.2 are ^ = 
1.46087 and = 1.93921; ^ = 1.48275 and 
^ = 2.41653 respectively (for T = 1). Hence 
Eqs. (I33i!l and (l33bll can be regarded as an 
approximate formula for constructing Arnold 
tongues. 

The faster is convergence of Feigenbaum 6 
exponent, the closer — and — are to 1. Hence 
one can see that there is another reason for 
investigating the behavior of the Feigenbaum 
exponents. 

We can see from Figs. [5] and [6|, that for the 
cases Q = 1.2 and Q = 0.8 there are three 
critical temperatures T„, T^-^ and T^^, at which 
in absence of external magnetic field transition 
between different phases takes place: 

1. T„ is the temperature corresponding to para- 
magnetic - anti- ferromagnetic phase (P-2AF1) 
transition (known as Neel temperature). The 
dependence on Q is shown in Fig. [TWa): 

2. is the one corresponding to anti-ferro- 
magnetic - four - periodic modulated phase 
(2AF1-2M2) transition [Fig. Hl^b)]; 
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3. Tc2, the temperature at which transition 
between four - periodic and eight - periodic 
modulated phases (2M2-2M3) occurs (the Q 
dependence of this temperature quahtatively 
is the same as of TcJ. 

Investigating the behavior of the Lyapunov 
exponent we found that critical temperatures 
and T^.^ exist not for all Q in the area 
1 < Q < 2 (the fact that Lyapunov exponent 
does not become zero indicates the absence of 
bifurcation, hence the second order phase tran- 
sition). After Q = 1.5681 there is no more crit- 
ical temperature like T^^ and after Q = 1.5115 
there is no more critical temperature like T^g. 

(a) 



10 - 



Q 

(b) 




Q 

Figure 12: The plot of (a) r„ (Neel temperature) 
versus Q for J = — 1, 7 = 3; (b) Tc^ (critical tem- 



perature corresponding to transition between anti- 
ferromagnetic and four periodic modulated phases 
in absence of external magnetic field) versus Q for 
J = -1,7 = 3. 

Figure [5] shows that there is some Tc (upper 
bound temperature), below which the system 
undergoes P-2AF1-P phase transition between 
paramagnetic and anti-ferromagnetic phases, 
which exists only in the area Q > 1 (e.g. for 
Q = 1.2, J = -1, 7 = 3 = 5.2606340). 

8 Conclusion 

In this paper we investigated Q-state Potts 
model in an external magnetic field on the Bethe 
lattice for non-integer Q. The model was ex- 
actly solved by means of recurrence relation 
technique and a one-dimensional rational map- 
ping was obtained. 

We pointed the relation between phase tran- 
sition of second order and bifurcation points. 
The phase structure of the model is investi- 
gated by constructing Arnold tongues. For 
given statistical model they correspond to the 
area of existence of different phases, i.e. to 
phase diagrams. We have obtained analogs of 
Arnold tongues with winding numbers w = 1/2 
and w = 2/4 in period doubling regime. We 
found the transitions P-2AF1, P-2AF1-2M2- 
2AF1, P-2AF1-2M2-2M3-2AF1 in the case of 
fixed H and P-2AF1-P at fixed T. The be- 
havior of the mapping describing the system is 
sensitive to the values of Q: case < Q < 1 in 
contrast to 1 < Q < 2. We also observed the 
Lyapunov exponent that describes the large 
variety of phase transitions that occur in the 
Potts model. We compared chaotic regions 
varying temperature and magnetic field and 
pointed that in the first case it was richer. 

We studied the convergence of the Feigen- 
baum a and 6 exponents for the above men- 
tioned rational mapping for < Q < 1 and 
1 < Q < 2 separately. It was shown that for 
a real statistical system these exponents coin- 
cide with the famous universal ones with high 
accuracy. 
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We constructed the Arnold tongue analogs 
with winding number w = 1/3 in the three 
periodic window and calculated the a and 5 
exponents for this area. The obtained values 
justified the universality hypothesis once more. 

We have found an approximate method for 
constructing Arnold tongues through Feigen- 
baum exponents, based on universality hypoth- 
esis. The main advantage of this kind of ap- 
proach is the simplicity of numerical calcula- 
tions in contrast to the rigorous treatment of 
the problem. 
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Table 1: Feigenbaum exponents for Potts-Bethe mapping [Eq. (^] for T = 1; J = —1; 7 = 3. 
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